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ABSTRACT 

An exponentially expanding Universe, possibly governed by a cosmological con- 
stant, forces gravitationally bound structures to become more and more isolated, even- 
tually becoming causally disconnected from each other and forming so-called "island 
universes" . This new scenario reformulates the question about which will be the largest 
structures that will remain gravitationally bound, together with requiring a systematic 
tool that can be used to recognize the limits and mass of these structures from obser- 
vational data, namely redshift surveys of galaxies. Here we present a method, based 
on the spherical collapse model and A^-body simulations, by which we can estimate 
the limits of bound structures as observed in redshift space. The method is based 
on a theoretical criterion presented in a previous paper that determines the mean 
density contrast that a spherical shell must have in order to be marginally bound 
to the massive structure within it. Understanding the kinematics of the system, we 
translated the real-space limiting conditions of this "critical" shell to redshift space, 
producing a projected velocity envelope that only depends on the density profile of 
the structure. From it we created a redshift-space version of the density contrast that 
we called "density estimator" , which can be calibrated from A-body simulations for a 
reasonable projected velocity envelope template, and used to estimate the limits and 
mass of a structure only from its redshift-space coordinates. 

Key words: methods: A-body simulations - large-scale structure of Universe - 
galaxies: clusters: general - galaxies: kinematics and dynamics 



1 INTRODUCTION 

The current paradigm of an exponentially expanding Uni- 
verse implies that large-scale, gravitationally bound struc- 
tures will eventually become causally disconnected from 
each other, forming island universes scattered inside a 



mostly empty Universe (e. g., Adams & Laughlin 1997 Chi- 
ueh fc He|2002l|Nagamine fc Loeb[20b3, |Busha et al. [2003 1 



In our previous paper (Diinner et al. 2006 hereafter "Pa- 



per I"), we presented a criterion to determine the limits 
of bound structures, defining superclusters as the biggest 
gravitationally bound structures that will be able to form. 
This criterion defined a critical density contrast over which 
a spherical shell will stay bound to a spherically distributed 
overdensity. As defined, this criterion can only be applied to 
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data given in real, three-dimensional space ("real space"). 
This is not the case of observational data, which comes from 
large galaxy surveys having two angular coordinates and a 
velocity coordinate. Using the Hubble law, one can estimate 
the real distance to an object from its recession velocity, 
but, given that we are interested in dense, relatively evolved 
structures with significant peculiar velocities, our estimation 
will be strongly affected by the velocity dispersion of the 
structure, fooling any attempt to apply a real-space-based 



method (Kaiser 19871 



Here we present a way to apply our theoretical criterion 
to redshift-space data, permitting its application to redshift 
surveys. The new criterion is based on the geometrical ap- 
pearance of the real-space criterion as seen in redshift space, 
and needs to be calibrated using statistics from A-body sim- 
ulations to account for velocity dispersions not considered 
in the theoretical model presented in Paper I. This method 
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represents an alternative to the caustic approach, first pro- 



posed by Regos & Geller ( 1989 1 and further developed by 
Diaferio & Gellerl il997\ and iDiaferiol (|1999|). The caustic 



method, based on the direct search for caustic curves which 
represent the redshift-space envelope of the bound structure 
within it, has been extensively used to study galaxy clus- 
ters (e. g., van Haarlem & van de Weygaert 1993 Geller, 



Diaferio, & Kurtz 1999, Reisenegger et al. 2000; Rines et al. 



2000, 2002, 2003: Biviano & Girardi,.2003^ Diaferio, Geller, 



& Rinos 2005, Rin es fc Diaferio|2006 1, constituting the most 
widely used method in the area. Among its main achieve- 
ments, it has been used to measure the mass profile and light 
to mass radio from galaxy clusters. 

Our method, even though it shares some of the basic 
elements of the caustic method, as will be discussed later, 
has the independent motivation of directly representing the 
spherical collapse density criterion for the critical shell in 
redshift space, giving a clear physical interpretation to its 
results. 

In Section [2] we discuss the effects of transforming the 
real-space data into redshift space as seen from A^'-body sim- 
ulations, introducing the critical projected velocity envelope, 
which is a theoretical construction produced by joint projec- 
tion of all shells, within a certain radius, intersecting the line 
of sight. We also study the implementation of a parameter- 
ized template for the density profile, for which we used the 



density profile proposed by Navarro, Prenk & White ( 1997 1 



In section |3] we propose a new method for applying 
the criterion in redshift space, which is presented in three 
alternative versions. For this we introduce the concept of 
"density estimator" , which replaces the previously used den- 
sity contrast in determining the threshold that defines the 
critical shell. We statistically determine the value of this es- 
timator, together with analyzing the main divergences from 
the spherical collapse theoretical model that introduce errors 
into the method. Each version of the method will require one 
or more density estimators, which will be presented later. 

In section|4] we test our method using our A'^-body sim- 
ulation, estimating the systematic error that is expected for 
radius and mass estimations for the gravitationally bound 
structure. 

Finally, in section [5] we present our conclusions, to- 
gether with a step by step recipe for applying one of the 
proposed methods. 



2 INGREDIENTS FOR A FITTING METHOD 

2.1 Redshift-Space Representation of the 
Spherical Collapse Model 

In Paper I, we showed that the spherical collapse model, 
when extended to the case of a universe dominated by a 
cosmological constant, can be used to set a criterion for the 
"critical" (marginally bound) shell of a mass concentration. 
The spherical shells are characterized by a single parameter 
named density contrast ils, where the s stands for shell. Its 
value, for the critical shell (cs), can be written as 



^cs — 



Pc 



= 2.36, 



(1) 



where pm is the mean mass density enclosed by the shell, 
and Pc = 3H'q/8tvG is the critical density of the Universe. 
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Figure 1. Three-dimensional effect of transforming to redshift 
space using the Hubble equivalence between velocity and distance. 
On the left, structure in real space. On the right, structure in red- 
shift space. In green particles bound to the structure and in blue 
unbound particles. For clarity, only particles inside the critical 
radius in real space were plotted. 



In the simulations, this criterion was shown to give 
an external limit to the extension of gravitationally bound 
structures, overestimating their mass by 39% on average. 
Nonetheless, the model gives in-falling velocity predictions 
which correctly follow the lower envelope of radial veloci- 
ties deep into the virialized core of bound structures (see 
Fig. 8 in Paper I). This, together with the one-to-one rela- 
tion between the shell's infall speed and its enclosed density 
contrast, makes it possible to extend the model to a redshift- 
space representation. 

To go from the real-space to the redshift-space repre- 
sentation, we need to replace the coordinate along the line 
of sight by the corresponding recession velocity. For sim- 
plicity, we will assume that the distance to the structures 
is much greater than their size, so we do not need to ac- 
count for angular effects and the replacement can be done 
directly in Cartesian space. An example of this transforma- 
tion is shown in Figure [l] where the structure on the right 
is the redshift representation of the structure on the left, as 
seen by an observer looking along the Z axis. In the case 
of a spherically symmetric structure, as described by the 
spherical collapse model, the projected velocity seen by an 
observer will be composed by the projection of the radial 
velocity of the spherical shell with respect to its center and 
the recession velocity of the whole structure (see Figure[2]A_). 

Let us consider a spherically symmetric, gravitationally 
bound structure. Each shell has its own radial speed, which 
only depends on its enclosed density contrast. As the spher- 
ical collapse model constrains the shells not to cross each 
other, we will see the innermost shells falling at greater 
speeds than outer ones. When moving from the center to 
greater radii, we will cross a shell that just stopped its ex- 
pansion, through slowly expanding ones, up to shells ex- 
panding with the Hubble flow. According to this model, the 
critical shell will be expanding with a speed (critical veloc- 
ity) of only 29% of the Hubble flow at the present time (see 
Paper I). 
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We are interested in identifying all shells lying within 
the critical shell, each of which will have a different projected 
shape in redshift space. Slowly expanding shells will appear 
as ellipsoids shrunk along the line of sight, while fast con- 
tracting ones will appear elongated along the line of sight, 
but flipped, so that apparently closer points will really be 
at the more distant side of the structure (Kaiser 119871. A 
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diagram explaining this effect is shown in Figure |2| where 
a fast contracting shell has a higher projected velocity than 
the outer expanding shell, protruding from the outer ellip- 
soid and producing the well-known "Finger of God" effect. 

We define the projected velocity envelope as the surface 
that encloses the redshift-space representation of all shells 
within a given radius R. This can be written as 



■ (rp) = max 

re[rp,fl] 



\Vr (r) 




(2) 



where is the radial velocity respect to the centre of 
the structure and rp is the projected radius on the sky 



i Reisenegger et aLl|2000 i. The spherical collapse model al- 



lows to predict Vi{r), and therefore v™^ (rp), directly from 
the density profile (Paper 1). Setting R equal to the critical 
radius of the structure (res), yields the critical projected ve- 
locity envelope ("critical envelope", for short), correspond- 
ing to the redshift-space envelope of all shells within the 
critical shell. The velocity envelope, when defined up to the 
turn-around radius, is equivalent to the caustics defined by 



Regos & Geller ( 1989 1 



If all the assumptions of the spherical collapse model 
were true, one would expect that all the bound particles 
should lie within this critical envelope. Moreover, if inner 
ellipsoids protrude from the ellipsoid defined by the critical 
shell, then one would expect to see contamination in the 
protruding regions from objects outside the critical shell. 

An effect that is not considered in the spherical collapse 
model are the velocity perturbations due to the interaction 
between nearby objects that will mix many particles into 
and out of the critical envelope, causing large systematic and 
random errors in the intended criterion for determining the 
limits of the bound structure in redshift space. To account 
for these systematic errors we used A''-body simulations as 
described in the next section. 



2.2 Velocity Envelope in A*'-Body Simulations 

In order to study the behaviour of structures in redshift 
space, we used numerical simulations, which permitted us 
to observe the redshift-space distribution of the structures, 
knowing at the same time which particles were bound and 
which were not. Our simulations are the same as described 



in Paper I, performed with the GADGET2 code (Springel 



2005), containing 128 dark matter particles inside a box 



of side length 100 Mpc, and considering a flat ACDM 
universe with Qa ~ 0.7. We took snapshots at the present 
time (a — 1) and in the far future (a — 100), assuming that 
in late epochs the structure evolution will decrease signifi- 
cantly so no major changes will be seen from then on (see 
also Busha et al.|2003 and Nagamine fc Loeb|2003l ). 

We selected the 11 largest structures for our study, with 
masses ranging from 1 X 10" Mq to 7 X lO" Mq. Although 
these structures might be rather small to represent our cur- 
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Figure 2. Redshift-space appearance of spherical shells. Panel 
(a) shows the position of shells in real space: the interior shell 
(grey) is contracting, the dashed shell is at rest, and the outer 
shell (black) is expanding at the critical velocity. The black dot 
on the contracting shell is falling faster than the expansion of 
the white dot on the critical shell, so its projected velocity is 
higher, appearing closer in redshift space. Panel (b) shows the 
corresponding ellipsoids in redshift space: the colouring has been 
kept from panel (a). The critical shell appears shrunk along the 
line of sight because its expansion speed is less that the Hubble 
flow at its radius. The interior shell instead appears elongated 
along the line of sight because its contraction speed is higher 
than the Hubble flow at its radius, but also flipped, such that the 
far side appears closer that the near side. The zero velocity shell 
appears as an horizontal line. The black long-dashed line shows 
the projected velocity envelope obtained from including all shells 
inside the critical shell. 



rent understanding of superclusters, many of them showed 
significant substructure, as expected from objects that are 
still evolving into a virialized state. The bound particles were 
identified using the state at a = 100 and then correspond- 
ingly tagged and followed to the present frame, repeating 
the procedure described in Paper I. 

In order to transform the simulated data to redshift 
space, we replaced the distance along one axis by the corre- 
sponding projected velocity. This can be done in any direc- 
tion, but for simplicity we did so only along the three main 
axes, giving a total of 33 data sets for statistical analysis. 

Figure |3] shows the critical envelope for two objects in 
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Obj #2 

Obj #3 

■■ Obj #e 

Obj tlO 



a 1 




Table 1. Mean values and standard deviations for several per- 
formance indicators for the critical envelope corresponding to the 
true density profile (see text). 



Figure 3. Critical projected velocity envelopes along the line of 
sight as expected in the spherical collapse model from the density 
profile for four well-behaved, representative simulated structures. 
Badly behaved structures, generally implying double cores or very 
significant substructure, usually present lower peak velocities and 
wider central regions. Radii are normalized in terms of the crit- 
ical shell's radius and velocities in terms of the Hubble flow at 
that radius. The masses satisfy Mahj#2 > -'^obj#3 > A'/obj#9 > 



(see Table |3|. 



obj#l 



our study with nearly spherical symmetry, with radii nor- 
malized in terms of the critical radius res and velocities re- 
ferred to the Hubble flow at res ■ These critical envelopes were 
obtained measuring the density contrast at all radii and ap- 
plying the procedure described in Paper I to find the radial 
velocity, and finally projecting it along the line of sight. At 
large radius, where there is no ellipsoid crossing, we see the 
ellipsoid corresponding to the critical shell, which has the 
same shape since its expansion velocity is the same fraction 
(29% in present time) of the Hubble flow (see Paper I). At 
smaller radius we observe the existence of ellipsoid cross- 
ing in redshift space, which is produced because the inner 
shells are contracting faster than the expansion velocity of 
the critical shell. Adjusting the object's centre, we observe 
an improvement in matching the resulting profiles. Specif- 
ically, a choice closer to the densest core of the structure 
increases the height of the peak of the velocity envelopes 
and improves the agreement between the shapes of curves 
corresponding to different objects. Even though the profiles 
do not match exactly at small radii, all of them share the 
same characteristic shape, showing higher velocities in less 
massive objects, meaning a higher concentration than more 
meissive structures, in qualitative agreement with the results 



of Navarro, Frenk & White ( 1997 1. These properties suggest 



that it may be possible to estimate the critical envelope by 
a general template that depends only on the bound mass of 
the structure. 

To study the accuracy of the prediction given by the 
critical envelope for the location of particles in redshift 
space, we counted the number of bound particles inside and 
outside the critical envelope, as well as the number of un- 
bound particles inside the critical envelope. Our results are 
summarized in Table [l] We used the same statistical indica- 
tors as in Paper I, so now we can compare the performance 





Redshift Space 


Real Space 


Indicator 


Mean 


Std. Dev 


Mean 


Std. Dev 


A 


23.96 


10.01 


28.2 


13.0 


B 


76.04 


10.01 


71.8 


13.0 


C 


30.59 


13.24 


0.26 


0.23 


D 


106.62 


17.56 


72.0 


13.1 



of the criteria in redshift and real spaces. The indicators, 
all expressed as percentages of the total number of particles 
mside the critical envelope, are: 

• A: particles inside the critical envelope that do not be- 
long to the cluster at a = 100. 

• B: particles inside the critical envelope that do belong 
to the cluster at a = 100. 

• C: particles outside the critical envelope that belong to 
the cluster at a = 100. 

• D: all particles that belong to the cluster at a = 100. 

Note that A + B = 100% and D = B -I- C. 

We observe that, compared to the real-space criterion 
(Paper I), the new redshift-space criterion produces better 
results for the indicators A and B, but the indicator C in- 
creases significantly (see Table [TJ , implying that it does not 
give an external limit to the location of particles in redshift 
space, in contrast with the real-space case. 

As pointed out before, compared to predictions done in 
real space, predictions done on redshift space are affected 
by a higher mixing of bound and unbound particles. This is 
due to the velocity dispersion, which makes many particles 
(inside the critical shell in real space) fall outside of the 
critical envelope in redshift space. The velocity dispersion 
is produced by local interactions between nearby objects as 
they decelerate driven by the central attractor, manifesting 
as local peculiar velocities in random directions. To check if 
the velocity dispersion was accompanied by an increase in 
the mean speed of the particles (as expected), we plotted the 
later with respect to their radii, comparing it to the absolute 
value of the theoretical radial velocity profile. In Figure [Ij 
we can see how the theoretical profile approximately marks 
the lower bound of absolute velocities for particles inside 
the critical radius, but outside the virialized core (around 
2.5 Mpc from the centre), so we can claim that the main 
contributor for particles escaping from the critical envelope 
is the overall increment in particle speeds due to interactions 
between them. 



2.3 Density Profile Template 

Until now, we have used the measured density profile to es- 
timate the critical envelope from our simulated structures. 
When confronted with observational data in redshift space, 
we will not have this information, so we need a general way 
to estimate the density profile of a structure based on a small 
number of parameters. We have seen in §2.2| that the criti- 
cal envelopes for different structures have a coherent shape 
(see Fig.|3|, but with scales correlated with the bound mass 



Redshift- space limits of bound structures 5 



1800 










1600 










1400 








- 


1200 


ml- 




1 


- 


1000 










800 
600 
400 






'f ' 


200 


r- - - ;- . -^Jy' 



















: 2 


4 


6 8 


10 1 



r [Mpc ] 

Figure 4. Absolute value of the velocities versus their radii to the 
centre for Obj. #9. The solid line shows the absolute theoretical 
velocity. The dashed line indicates the critical radius. The turn- 
around radius is clearly observed where the solid line touches zero, 
accompanied, as expected, by a minimum of absolute velocities. 
For clarity we plotted a 70% random sample of all particles. 



of the structure. This is in agreement to what was pro- 



posed by Navarro, Frenk & White (19971, who presented 



the now well-known NFW profile, which is a generalized 
density profile for virialized clusters, obtained empirically 
from many simulations, fn particular, we used the results 
from a later publication ( Eke, Navarro, &: Steinmetz|2001 1, 
which corrected several details from the original work. In 
general words, they postulate that a cluster's density pro- 
file can be estimated using as a single parameter, the radius 
r A at which a characteristic ( "virial" ) density contrast A is 
reached. The parameter A depends on the cosmology and 
can be obtained in a fiat universe as A = n8^l°J\ which m 
our case gives A = 103.5. This characteristic radius is much 
smaller than our critical radius, falling near the edge of the 
virialized zone where clusters are less affected by external 
structures and look more homogeneous, making it easier to 
estimate using redshift data. The NFW profile was formu- 
lated to give good results for radii ranging from O.Irs to lOr^, 
where — va/ca and ca is the concentration parameter, 
which can be obtained as a function of the mass contained 
within rA, and whose value is around ~ 8 for the masses of 
the structures studied by us. Considering that the critical 
radius is of order 3 — 5rA, the NFW profile is more accurate 
for the inner part of the structure. This region is the most 
interesting one for our purposes, since it is where ellipsoid 
crossing takes place and the shape of the critical envelope 
is directly dependent on the density profile. On the other 
hand, we caution that large bound structures may generally 
not be dominated by clusters as relaxed and spherical as 
those for which the NFW profile was derived. 

The NFW density profile was fitted to every object in 
our simulation using the actual density in real space. 

Figure |5] shows the density profiles for all our simulated 
structures, compared to the NFW template profile. Scales 
have been normalized by their corresponding NFW scale, in 
order to let us compare all profiles to a single NFW tem- 
plate. The dark and light grey vertical areas show where the 
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Figure 5. Density profiles of all structures (solid), compared to 
the NFW template (dashed). The densities and radii have been 
normalized by their corresponding scale in the NFW template to 
put them on the same plot. The dark and light grey vertical areas 
denote placement of critical radii measured from the simulated 
data and from the fitted NFW template respectively (the dark 
grey area actually extends under the light grey area up to its left 
end, but is hard to see in the plot). The horizontal area denotes 
the range for the asymptotic values for the density, QrnPc- 



critical radii obtained from the simulated data and from the 
NFW profile, respectively, are located. The horizontal area 
shows the asymptotic values expected as the density reaches 
the mean density of the Universe. Clearly, the true density 
profiles depart from the template at radii somewhat smaller 
than the critical radius, as they approach the mean density 
of the Universe. This early departure from the model implies 
that estimation of the critical radius done directly from the 
NFW template will be biased to lower values. Below, we will 
present ways to deal with this bias. 

As seen in Figures|6]and[7[ the radial velocities and crit- 
ical envelopes predicted using a NFW density profile show 
different levels of agreement with the true velocities depend- 
ing on the range of radii considered. For small radii, it yields 
a higher density than observed, predicting higher infall ve- 
locities. Given the resolution of our simulation, which was 
intended to search for large scale structure rather than repli- 
cating the behaviour of the virialized cores of structures, we 
believe that our simulated data is unable to produce ac- 
curate densities at such radii. For intermediate radii, the 
prediction is very good, as expected, since the NFW pro- 
file was fitted at a measured rA in this range. Finally, for 
larger radii, closer to res, the generalized profile gives better 
or worse results depending on the absence or presence of sig- 
nificant substructure, respectively. In general, the presence 
of substructure increases the density at higher radii, so the 
NFW profile underestimates it, consequently giving a low 
estimate for the critical radius. 

Tables [2] and |3] compare the predicted critical radius 
and enclosed mass for every object in the analysis. We ob- 
serve that the critical radius according to the NFW profile 
is 89.8% of the measured critical radius. Thus, considering 
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Figure 6. Radial velocities compared to the theoretical velocity profile obtained from the true density profile (black solid line) and from 
the NFW density profile fitted to the structure (red dashed line). Green dots indicate bound objects and blue dots the opposite. On the 
left object #1, and on the right object #9. Sampling 50%. 




Figure 7. Projected velocities versus normalized projected radius. In green, particles that will eventually fall in the object. In blue, 
particles that will escape. The curves show different critical envelope estimations: The black solid line shows the true velocity envelope. 
The red dashed line is the NFW envelope adjusted using the observed t-a value. The cyan solid line is the NFW envelope fitted using 
the NFW-core method (fitting methods are explained in ^3.2[ |. The magenta solid line is the NFW envelope fitted using the NFW-cs 
method. The brown solid line show the combined envelope fitted using the combined method. On the left object #1, and on the right 
object #9. SampUng 50%. 



that the spherical collapse criterion gives an external limit 
for the critical radius, estimations done with the NFW pro- 
file will underestimate the size of the structure as defined by 
the spherical-collapse criterion. 

Concerning the mass of the structures, we find that the 
mean true bound mass (Mbound) is 82.7% of the mass en- 
closed by the NFW profile critical radius in real space. This 
should be compared to the same relation for the true criti- 
cal radius, where the true bound mass is 71.7% of the mass 
enclosed by it (see Paper I for details). 



3 FITTING METHOD DEFINITION 
3.1 Density Estimator 

The great advantage of using the NFW profile is that the 
resulting critical envelope depends exclusively on the bound 
mass of the studied structure. In this way, the projected ve- 
locity envelope can be scaled (changing either or Ma) 
until the contained mass and volume satisfy some condi- 
tion equivalent to the critical density contrast, but in red- 
shift space. For this purpose, we will define an observable to 
characterize the redshift-space density contrast inside the 
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Table 2. Critical radii considering different models in real space. 
Columns: (1) Object index. (2) r^: virialization radius needed to 
fit the NFW density profile. (3) res: critical radius applying to 
our criterion {Qcs = 2.36) to the true density profile. (4) r^^*: 
critical radius applying our criterion to the fitted NFW profile 
(QNFW = 2.36). (5) Ratio between r^/W ^nd res. 



Obj.# 


r A [Mpc] 


res [Mpc] 


reN/W[Afpc] 


„NPW /„ 
'cs /'cs 


1 


1.81 


8.28 


8.21 


0.99 


2 


1.51 


7.84 


6.79 


0.87 


3 


1.48 


7.54 


6.65 


0.88 


4 


1.56 


7.92 


7.02 


0.89 


5 


1.30 


7.38 


5.80 


0.79 


6 


1.21 


7.93 


5.38 


0.68 


7 


1.48 


7.63 


6.64 


0.87 


8 


1.30 


6.29 


5.80 


0.92 


9 


1.31 


6.01 


5.85 


0.97 


10 


1.27 


5.67 


5.66 


1.00 


11 


1.23 


5.37 


5.49 


1.02 



Table 3. Enclosed masses considering different models in real 
space, contrasted to the effectively bound mass. Columns: (1) 
Object index. (2) M^ound- Bound mass at a = 100. (3) Mes: 
Mass enclosed by the critical radius res. (4) Mass ratio between 
bounded mass and Mes (5) M^^^^: Mass enclosed by the critical 
radius obtained using the NFW density profile r^^^. (6) Mass 
ratio for M^^W, 



Obj.# 


-^-^bound 


Afes 


«bound 


A/e^FW 


A/bound 


1 


6.34 


7.64 


0.83 


7.60 


0.83 


2 


5.29 


6.47 


0.82 


5.96 


0.89 


3 


4.71 


5.76 


0.82 


5.26 


0.90 


4 


4.14 


6.68 


0.62 


5.99 


0.69 


5 


4.07 


5.41 


0.75 


4.31 


0.94 


6 


3.23 


6.70 


0.48 


3.01 


1.07 


7 


2.95 


5.98 


0.49 


4.53 


0.65 


8 


2.44 


3.34 


0.73 


3.18 


0.77 


9 


2.40 


2.91 


0.83 


2.79 


0.86 


10 


2.02 


2.45 


0.82 


2.44 


0.82 


11 


1.46 


2.08 


0.70 


2.18 


0.67 


All masses in units of 10^* 









surface given by a certain projected velocity envelope. The 
redshift-space "density estimator" is defined as 



where Mr is the mass inside the projected velocity envelope 
defined by the radius R and Vr — iivR^ /3 is the real-space 
volume enclosed by the same radius. This definition allows 
us to obtain density estimators for any projected velocity 
envelope, including the critical envelope. In practice we will 
need to define density estimators for only two kinds of pro- 
jected velocity envelopes: the critical envelope and the "core 
envelope" , which is the projected velocity envelope obtained 
using the NFW density profile for radii within r a , necessary 
to estimate the NFW density profile in redshift space. These 
density estimators can be computed for all the structures in 
our sample, and then averaged to produce a calibrated es- 
timator that can let us infer the correct critical envelope 
without need of real-space data. 



The aim is to have a general template for the desired 
projected velocity envelope, which can be associated to a 
specific density estimator, specially calibrated for it. When 
fitting a template to data from a redshift survey, it can be 
scaled until its density estimator takes the desired value, 
yielding an estimation for the radius in real space that char- 
acterises it. 

To understand better the properties of the density esti- 
mator, we calculated it using the "true critical envelope" (in- 
fered from the true density obtained from the simulations) 
from all our 11 objects from 3 points of view. Results are 
shown in Table |4] The mean value for the density estimator 
is f2es = 1.60. This value is significantly lower than the cri- 
terion Qcs = 2.36 for real space, implying that the selected 
volume is less populated or, in other words, many more par- 
ticles escaped the theoretical prediction, also reflected in the 
greater value for C. This leads to the conclusion that per- 
turbations away from the spherical collapse model have a 
stronger manifestation in particles velocities than in their 
positions, decreasing the reliability of predictions in redshift 
space. 

3.2 Critical Envelope Recipes and Fitting Method 

We would like to replace the true density proflles (unknown 
in actual observed structures) by a simple parametrisation, 
based on the NFW density proflle. In order to test its accu- 
racy, we compare the previously defined performance indi- 
cators for the following three projected velocity envelopes: 

• True envelope: The projected velocity envelope is de- 
termined directly from the true density profile. 

• NFW envelope: The virialization radius ta is deter- 
mined using the appropriate over-density criterion in real 
space, yielding an NFW density profile, which is used to gen- 
erate a projected velocity envelope defined for shells within 
the critical radius according to the NFW density profile. 

• Combined envelope: rA and res are found using the ap- 
propriate over-density criteria in real space. We generate a 
projected velocity envelope using the ellipsoid correspond- 
ing to the critical shell, adding the inner part of the NFW 
envelope where it protrudes from the ellipsoid. 

Table [5] gives a comparison of the statistical indica- 
tors for all the three critical envelopes in consideration. We 
clearly see how the combined envelope shows almost the 
same statistical properties as the true envelope. In contrast 
the NFW envelope shows a greater amount of particles that 
escaped the profile (Indicator C). 

Until now, the projected velocity envelopes have been 
determined from information in real space that is not avail- 
able in actual observed structures. To fit the velocity en- 
velopes directly in redshift space, we will need to use density 
estimators in order to scale the templates up to the desired 
size. Here we propose three ways to do this. 

The first way is to directly try to estimate the NFW 
density profile for the structure, for which one needs to es- 
timate the virialization radius rA- In order to estimate ta 
from the redshift space data, we defined another template 
based on the NFW density profile: 

• Core envelope: equivalent to the NFW envelope, but 
now considering only shells inside the virialization radius. 
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Table 4. Density estimators for all projections. Rows: (1 - 3) Projections in the x, y and z axes. (4) Mean values for each object. 
Columns: All 11 objects. 



Object 


#1 


#2 


#3 


#4 


#5 


#6 


#7 


#8 


#9 


#10 


#11 


X axis 


1.92 


1.48 


1.65 


2.16 


1.23 


1.30 


1.43 


1.84 


1.77 


1.68 


1.68 


y axis 


1.54 


1.81 


1.63 


1.62 


1.34 


1.20 


1.37 


1.50 


1.77 


1.81 


1.67 


2 axis 


1.55 


1.58 


1.76 


1.45 


1.40 


1.12 


1.53 


1.49 


1.84 


1.83 


1.87 


Mean 


1.67 


1.62 


1.68 


1.74 


1.32 


1.20 


1.44 


1.61 


1.79 


1.77 


1.74 



Table 5. Mean values and standard deviations for several performance indicators for three critical envelopes obtained from real-space 
measurements. Rows: (1) True: Critical envelope calculated using the true density profile. (2) NFW: Uses the NFW density profile 
according to r/^ in real space. (3) Comb: Uses the NFW profile for the innermost radii and the maximum ellipse given by the critical 
radius. 



Envelope 


A 




B 




C 




D 




True 


23.96 


10.01 


76.04 


10.01 


30.59 


13.24 


106.62 


17.33 


NFW 


19.72 


7.35 


80.28 


7.35 


37.17 


21.67 


117.42 


24.23 


Comb 


23.88 


9.98 


76.12 


9.98 


30.80 


14.09 


106.93 


17.55 



Then we can construct a density estimator fi^^J^ that can 
be used to fit the core envelope in redshift space. Once we 
have estimated ta, we can use it for estimating r^^^, de- 
fined as the radius where the density contrast in real space 
^NFw ^ 2.36. We labeled this recipe as "NFW-core". 

The second way is to directly fit the NFW envelope, 
which is defined up to r^^^ , using a density estimator pre- 
viously calibrated for this template (fi^^^). The fitted en- 
velope will directly estimate the critical radius r^^^ . We 
labeled this recipe as "NFW-cs" . 

The last method is to use the combined envelope, for 
which we need to estimate two independent parameters, ta 
and Tcs- For the first we can fit the core envelope, following 
the same procedure as in the NFW-core recipe. After esti- 
mating TA, we can construct the combined envelope leaving 
as a single free parameter the critical radius res, which de- 
fines the outer ellipsoid. Finally we can use another density 
estimator, previously calibrated for the combined en- 

velope, to fit this last parameter. We labeled this recipe as 
"comb" . 

Having estimated the critical radius using one of these 
methods, we use equation ([T]) to estimate the enclosed mass. 
Then the true bound mass can be estimated using the statis- 
tical relation between the mass inside the critical shell and 
the bound mass, as explained in Paper I. 



3.3 Centre Determination 

A very important step in our method for identifying the 
structure limits is to first have a good estimate for its cen- 
tre, where to anchor the projected velocity envelope. The 
intuitive way to do this is to look for the region with the 
highest density of particles (galaxies) in redshift space. In 
other words, we need to find the centre which maximises 
the size of the projected velocity envelope for a given den- 
sity estimator. We performed this maximisation using the 
core envelope, focusing on finding the densest virialized cen- 
tre for our structure. As density estimator we used f2^rc* 
calibrated for the centre found in real space, using the pro- 
cedure described in the next subsection. 



We found that the new centre appeared to be displaced 
from the original centre on average a 13.5% of va (represent- 
ing the virialized core), but the errors produced from using 
this new centre were marginal compared to other systematic 
errors, being around one order magnitude smaller. 



3.4 Density Estimator Calibration 

In order to use the general projected velocity envelope to 
estimate the critical radius of gravitationally bound struc- 
tures, we first need to obtain a density estimator specially 
calibrated for that particular type of projected velocity en- 
velope. The best way of doing this is to use a statistically 
significant number of simulated structures in the scale range 
we are interested in, and calculate the corresponding den- 
sity estimator using the actual projected velocity envelope 
obtained from real-space measurements. 

Here, we are going to use the mean of the density es- 
timators found in all 33 projections as our calibrated den- 
sity estimator, and the standard deviation will be used as 
the expected error. The density estimators needed by our 
projected velocity envelope recipes are rjJ^^J^ for the core 
envelope, fi^/^ for the NFW envelope and 0.1°'^^ for the 
combined envelope. 

The mean density estimators, together with their stan- 
dard deviations are shown in Table [6] For completeness, and 
as a comparison point, we also included the mean density 
estimator for the true envelope, f^cs"'^, which can also be 
found for every single object and point of view in Table |4] 

We observe that the True and Combined density esti- 
mators are very similar, which tells us that the NFW profile 
accurately reproduces the central part of the true velocity 
envelope. On the other hand, the NFW envelope has a much 
greater density estimator, which is the result of a smaller 
critical radius (r^^^). The density estimator associated to 
the core envelope is actually very close to its equivalent value 
in real space (A = 103.5), which is good considering that 
the structures' core are ruled by virialization, far from the 
spherical collapse constraints we are assuming. 
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Table 6. Mean values and standard deviations for inner density 
estimators for several projected velocity envelopes with parame- 
ters determined in real space. See row definitions in the text. The 
first three estimators correspond to envelopes defined up to the 
critical radius, while the forth is defined up to the much smaller 
virialization radius r/\. 



Envelope 


Estimator 


Mean 


Std. Dev. 


True 




1.59 


0.25 


NFW 




2.06 


0.35 


Comb 




1.59 


0.26 


Core 


qNPW 


109.64 


14.56 



3.5 Main Sources of Error 

We can distinguish two kinds of error sources when first 
looking at Table [4] the first is associated to big variabil- 
ity in the estimator depending on the axis being viewed, 
and the other is manifested as an overall under or overes- 
timation of the density estimator for all views of the same 
object. For the first, the axis dependence tells us that there 
is a strong influence of the cluster's particular geometry or 
anisotropies outside of the spherical distribution. A good 
example is Obj.# 4, where we find a big deviation in obser- 
vations done along the X axis, showing an overestimation 
of the inner density with respect to observations done along 
the other axes. After analyzing the statistics, and plotting 
the radial velocities from all axes (see Figure|8|, we conclude 
that this difference is produced by an infiow of particles from 
an external structure placed right on the line of sight, and 
whose virialized velocities mix in redshift space with par- 
ticles from the studied cluster. The effect is an oversized 
parameter A and an undersized parameter C. 

These observations appeared repeated in every object 
with big standard deviation of the density estimator among 
the three axes, but did not appear related to objects with 
highly non-spherical cores or double cores. 

Regarding the second type of error, we conclude that it 
is related to highly non-spherical cores or anisotropies inside 
the critical radius. That is the case of Obj.#5 shown in Fig- 
ure [oj which has a double cor^ We observe that the double 
core does not produce a big deviation between estimations 
from different axes, but produces a significant underestima- 
tion of the inner density, which at the end will cause that 
this kind of objects will be bigger than predicted. 



4 FITTING METHOD TESTING 

Before being able to apply the method to real observations, 
we tested it in our simulations, in order to get statistical 
information that can be used to interpret later results. 

4.1 Performance Statistics 

We fitted the three kinds of critical envelopes to our data 
in redshift space using our previously calibrated density es- 
timators. Considering that we have only 11 objects in our 

^ In the analysis we used as centre the centre of the bigger core. 



study, we decided to use the same data to do the fittings, 
but we rotated the redshift projections in 45° to at least 
use different points of view at this stage. To fit the profiles, 
we just scaled them until the density estimator condition 
was reached. As centre for the profiles, we used our redshift- 
determined centres. We realized three fittings: 

• Combined method (comb): Uses the density estimator 
f^cOTc^ to fit the core envelope and calculate ta- Then fits 
an ellipsoid (that corresponds to the critical shell) adjusting 
its radius until the combined envelope produce the desired 
density estimator, f2J:°™'', directly estimating res- 

• NFW-based method (NFW-core): Uses the density esti- 
mator to fit the core envelope and calculate ta- Then 
uses the NFW density profile to calculate the critical radius 

NFW 
' cs 

• NFW-based method 2 (NFW-cs): Uses the density esti- 
mator f2^^^ to fit the whole NFW envelope to the structure 
in redshift space, directly estimating rj^^^. 

For the radius and mass estimations, we need to con- 
sider the systematic biases discussed in subsection |2.3| for 
the critical radius obtained from the NFW density profile 
r^sF^ and for the bound mass estimated from the mass 
enclosed by the critical radius. In particular we use that 
res = r^/W/0.898 and Afbound = 0.717Afcs, where M,, is 
the mass derived from equation Q. 

There are two main statistics we are interested in: The 
first denotes mean systematic deviations of the estimated 
parameter (radius or mass) from the true valu^ For this 
we define the parameter ai as the mean of the ratio between 
the estimated and the true value of the physical property i. 
The existence of this kind of mean biases is mostly due to 
statistical deviations, and should disappear as we increase 
the number of structures studied for the whole calibration 
and testing of this system. 

In second place we are interested in quantifying errors 
due to morphological characteristics that depart individual 
structures from the mean. A way to do this is calculating 
the standard deviation of measurements done from different 
angles, normalised by the true value, which has the advan- 
tage of discounting the mean bias produced by systematic 
errors. We labelled this as <Jx,i, where i denotes again the 
physical property we are referring to. The mean parame- 
ter, averaged over all structures, ax, will give us the error 
we expect to find in estimations done without taking into 
consideration the particular morphological properties of the 
studied structure. 

The mean results from this analysis can be found in Ta- 
ble [7] There we can observe that the systematic biases are 
all very small, reaching in the worst cases 1.5% for ta esti- 
mations, 0.7% for estimations and 8.6% for mass estima- 
tions. This result should be taken with some care, because 
we need to consider that there is an evident correlation be- 
tween the source and the test data for our system, but this 
situation can be improved very easily by growing our sample 
with new simulated structures. 

For the standard deviation, we observe that errors due 
to morphological anisotropies out of the spherical symmetry 



^ For the mass, the true value is the bound mass in a = 100. 
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Figure 9. Projected velocities versus normalized projected radius for Obj#5. This object presents a double core that will eventually 
merge into a single object. We observe that this feature does not produce a significant effect when observing it from a different point of 
view. Sampling 50%. 



Table 7. Systematic bias errors and mean standard deviation. 
Columns: (1): Fitting method. (2), (4) and (6): Mean normalised 
bias a in the measured , res and M. (3), (5) and (7): standard 
deviation , which represents the deviation of estimations along 
different axes with respect to the mean of the estimations for the 
same object. 



Method 


"A 












Comb 
NFWl 
NFW2 


1.013 
1.013 
1.015 


0.080 
0.080 
0.066 


1.000 
1.007 
1.003 


0.068 
0.083 
0.068 


1.045 
1.086 
1.055 


0.206 
0.253 
0.203 



reach in the worst cases 8.0% for ta estimations, 8.3% for 
res estimations and 25% for mass estimations. 



In Figures |10| |11| and |12[ we present the results for 
all methods. It is clear that the method is axis-dependent, 
showing big deviations in structures with strange morpholo- 
gies. These deviations are in agreement with the variability 
observed in the density estimators shown in Table|4] suggest- 
ing that much can be done identifying types of structures 
and calibrating custom density estimators for them. 



As a last conclusion, we observe that the best results are 
obtained with the combined and NFW-cs methods, both of 
which have in common that they are fitted to res instead of 
rA- 
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Figure 10. Estimated virial radii versus measured virial radii for all profile recipes and points of view. On the left, results from the 
combined method. In the centre, results from NFW-core method. On the right, results from NFW-cs method. Circles indicate X-axis 
projections, triangles indicate V-axis projections and squares indicate Z-axis projections. 



Figure 11. Estimated critical radii versus measured critical radii for all profile recipes and points of view. On the left, results from the 
combined method. In the centre, results from NFW-core method. On the right, results from NFW-cs method. Circles indicate Jf-axis 
projections, triangles indicate Y-axis projections and squares indicate Z-axis projections. Results have been debiased considering that 
J'cs.NFW = 0.898rcs 
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Figure 12. Estimated bound masses versus measured bound masses (in a = 1). On the left, results from the combined method. In 
the centre, results from NFW-core method. On the right, results from NFW-cs method. Circles indicate X-axis projections, triangles 
indicate Y-axis projections and squares indicate Z-axis projections. Results have been debiased considering that r^^^ = O.SQSrcs for 
the NFW-based methods and Mbound = 0.717Afcs for all. 



4.2 Error Estimation 

The most direct way to estimate the expected error from 
estimating the size and mass of a bound structure using this 
method is to calculate the actual deviations from the true 
value over all objects in our study. For the critical radius 
and the bound mass we used the following equation: 



Table 8. Expected errors in critical radius and bound mass es- 
timations. Columns: (1): Fitting method. (2): Percentage error 
in critical radius estimation. (3): Percentage error in bound mass 
estimation. 



Method 


Error,- % 


Error M% 


Comb 


8.5 


31.3 


NFWl 


12.2 


44.1 


NFW2 


8.1 


30.9 



where x is the parameter on which we want to determine 
the error. 



The expected errors calculated in this way are shown 
in Table [S] We emphasise that these errors are the ones 
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we should expect after a "brute-force" application of the 
method, having no consideration of the observable morphol- 
ogy of the structure. Due to the statistical nature of the 
fitting method, one would expect significant improvement 
in the expected errors after a more careful analysis of indi- 
vidual cases, where different "kinds" of structures, including 
structures with evident double cores or highly contaminated 
by substructure, were considered. This analysis is beyond 
the scope of this work, but can be included in future works. 



5 CONCLUSIONS 

We have presented a new method to estimate the mass 
and radius of gravitationally bound structures based solely 
on redshift information present in redshift survey catalogs. 
The method is based on the spherical collapse model (Pa- 
per I) and on the important observation that the theoretical 
critical envelope correctly follows the true critical envelope 
deep inside the virialized centre of the simulated structures. 



We used the NFW density profile ( Navarro, Frenk & White 
to generate a template for the critical envelope, which 
was calibrated using A^-body simulations. 

The extension of the method to redshift space gave birth 
to a new set of criteria that we called "density estimators" . 
These were defined as the expected redshift-space density in- 
side the projected velocity envelope. Their calculation was 
completely empirical, depending on statistical analysis of 
simulated structures due to the considerable velocity dis- 
persion found in them. We observed that the main cause of 
velocity dispersion was the gain of peculiar velocities due 
to local interaction between in-falling objects and substruc- 
ture. The study of substructure and of the morphology of the 
studied structure can be of great help to improve the error 
bars, since the density estimators show great dependence on 
these properties. From this, we conclude that numerical sim- 
ulations can be used to emulate the particular properties of 
the studied structures, producing "custom-made" methods 
to obtain the best results. 

In contrast to the caustic method, our method forces a 
limit to the radius of the structure, defined as the maximum 
radius at which one should expect to find bound objects. 
The shape of the velocity envelope though is less flexible 
than the caustics shape, suggesting that a combination of 
both methods can give better results, finding the envelope 
at lower radii using the caustics where they are more defined, 
and using the fixed envelope at higher radii to set a limit to 
the structure. 

The more reliable methods from the study were the 
"combined" and the "NFW-cs" methods. The procedure to 
apply the combined method to a redshift data set is then 
the following: 

• Calibrate data so that every element is accurately re- 
lated to a mass, so that the whole data set serves as a 
redshift-space mass field. 

• Identify the center of the structure by maximising the 
number of particles inside the core envelope for a given den- 
sity estimator. 

• Fit the NFW profile to the central region (core) by ad- 
justing ta until the measured density estimator is equal to 
the mean density estimator for the core, given in Table ID 



• Add the ellipsoid given by res to the critical envelope, 
and adjust res to the value fcs, at which the density esti- 
mator inside the total critical envelope is equal to the mean 
density estimator for the combined method, given in Table 

m 

• Given fcs, calculate the bound mass using equation [T] 

• The fractional errors for fcs and Mbound are given in 
Table H 

In our next paper we will present the application of the 
method presented here to estimate the size and bound mass 
of the Shapley supercluster ( Proust et al.|[2006 K 
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